Effect of annealing on the specific heat of Ba(Fei_ x Co x ) 2 As2 
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We report on the effect of annealing on the temperature and field dependencies of the low tem- 
perature specific heat of the electron-doped Ba(Fei_ :E Co a ;)2As2 for under-(x = 0.045), optimal- (x 
= 0.08) and over-doped (x = 0.105 and 0.14) regimes. We observed that annealing significantly 
improves some superconducting characteristics in Ba(Fei_ I Co a; )2As2. It considerably increases T c , 
decreases 70 in the superconducting state and suppresses the Schottky-like contribution at very 
low temperatures. The improved sample quality allows for a better identification of the supercon- 
ducting gap structure of these materials. We examine the effects of doping and annealing within a 
self-consistent framework for an extended s-wave pairing scenario. At optimal doping our data indi- 
cates the sample is fully gapped, while for both under and overdoped samples significant low-energy 
excitations possibly consistent with a nodal structure remain. The difference of sample quality offers 
a natural explanation for the variation in low temperature power laws observed by many techniques. 

PACS numbers: 74.20.Rp, 74.70.Dd, 74.62.Dh, 65.40.Ba 



I. INTRODUCTION 

It is not surprising that two years after the initial dis- 
covery of superconductivity (SC) in Fe-based materials^ 2 - 
there continues to be much debate as to the structure of 
the superconducting gap and the origin of the pairing 
mechanism. Even in the cuprates these questions are 
still not fully resolved, however significant progress was 
achieved with improved crystal quality. While the high- 
est T c is found in the 1111 family of compounds^, more 
progress on the gap structure has been made on the 122 
family of compounds^, as large single crystals can be 
synthesized. While much has already been learned on 
these crystals^, significant amounts of impurities are also 
clearly present as evidenced by the residual linear term 
in specific heal£~— , /xSR— , optical conductivit y 12 ' 13 , the 
broadening of the NMR lineshape with doping^ and the 
gap inhomogeneity observed by ST M 15 ' 16 . Some inho- 
mogencity may be unavoidable as alloying is currently 
required to achieve superconductivity at the highest tem- 
peratures. However, the availability of cleaner crystals 
could expose the intrinsic power law behaviors necessary 
to identify the superconducting gap structure and per- 
mit better phase sensitive experiments to be performed. 
In addition, the evolution in low temperature properties 
between clean and dirty samples can itself be a useful 
probe of the SC order parameter (see RefJl7l and ref- 
erences therein). For instance, the penetration depth 
of samples with varying amounts of heavy ion irradia- 
tion was recently studied^. Annealing is another typical 
route in the opposite direction: it improves crystal qual- 
ity. Indeed, for doped SrFc2As2 it was recently shown 
that annealing can significantly reduce scattering in the 
normal state, as well as enhance T ^ 19 ' 20 . 

Fe-based superconductors are widely believed to have a 



so-called s± gap structure, belonging to a fully symmet- 
ric A\ g representation but with a sign change of the order 
parameter between the electron and the hole Fermi sur- 
face sheets 2 ^ - — . While the gap on the hole Fermi surface 
(FS) sheets is believed to be nearly isotropic, there are 
indications from both theory 2 ^— and experiment* 2 ^— 
that on the electron sheet a strongly anisotropic, or 
even nodal, gap is produced by the competition of 
the magnetically-assisted interband scattering and the 
Coulomb interaction within each band (electron or hole) . 
The anisotropy is expected to be more pronounced on 
the overdoped side. 

Impurities are crucial for determining the precise low- 
energy density of states and, hence, the nature of the 
superconducting order parameter. First, the impurity 
scattering within each band tends to make the gap on 
the corresponding FS sheet more isotropic, and may re- 
sult in lifting the nodes and removal of the low-energy 
excitations In contrast, the nonmagnetic scatter- 

ing by the same impurities between bands with an op- 
posite sign of the superconducting order is pairbreaking, 
and increases the number of unpaired quasiparticles in 
the gap. For impurity concentrations exceeding a crit- 
ical value (which depends on the strength of the indi- 
vidual scatterers and the relative strength of the inter- 
vs intra-band scattering) the system possesses low en- 
ergy excitations. Therefore a study of these excitations 
and their correlations with impurity strength and con- 
centration addresses the nature of superconductivity in 
pnictides. 

One way to address these states is by studying the 
low-temperature heat capacity. The measurements of 
the temperature dependence of the specific heat and its 
magnetic response in the superconducting state give im- 
portant information about the symmetry of the order 
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FIG. 1: (Color online) The temperature dependence of the heat capacity of Ba(Fei_a; 00^)2 As2 before (red circles) and after 
(blue squares) annealing for: (a) x = 0.045, (b) x = 0.08, (c) x = 0.105 and (d) x = 0.14. Arrows mark the superconducting 
transitions. Inset: low-temperature part of the electronic specific heat of Ba(Fei_ I Co I )2As2. 



parameter— ~—. For near optimally Co doped BaFc2As2 
such an analysis suggest the presence of an anisotropic 
gap structure which may or may not have nodes (see 
ReflHHHil. As have been shown in the latter pa- 
pers, specific heat displays significant residual density of 
states as a consequence of impurities of some form (see 
also Refill). Moreover, the heat capacity at very low tem- 
peratures is very sensitive on the crystals quality and very 
often dominated by a Schottky-like contribution^. This 
can make the gap structure analysis quite complex. 

In this paper, we present results of our specific 
heat studies of the under-, optimally- and over-doped 
Ba(Fci_ x Coa;)2As2 as a function of annealing. Some 
preliminary results at optimal doping were presented in 
Ref® We show this process improves significantly su- 
perconducting characteristics in the Co-doped BaFc2As2 
(Bal22) system. It increases T c and results in a pro- 
nounced decrease of the residual linear term of the spe- 
cific heat. Moreover, after annealing, we do not observe 
a low-temperature Schottky anomaly. Subsequently, we 
examine the low temperature and magnetic field behav- 
ior of the electronic contribution to the specific heat and 
compare with self consistent calculations based on an ex- 



tended s-wave pairing mechanism with impurities. The 
suppression of low energy excitations as a function of an- 
nealing in the optimally doped sample suggests a fully 
gapped superconducting state. In contrast, the evolu- 
tion with annealing at under and over-doping indicate 
that the superconducting gap structure contains nodes 
at these doping levels. This implies a change in the gap 
structure with doping. 



II. EXPERIMENTAL DETAILS 

Single crystals were grown out of FeAs flux with a typ- 
ical size of about 2x1.5x0.2 mm 3 (see Ref H). They crys- 
tallize in well- formed plates with the [001] direction per- 
pendicular to the plane of the crystals. The doping levels 
were determined by microprobe analysis. After synthesis 
all samples have been characterized and studied, then we 
annealed the samples for two weeks at 800 °C in vacuum 
and repeated the measurements on the same pieces. The 
heat capacity was measured down to 400 mK and in mag- 
netic fields up to 9 T using a thermal relaxation method 
implemented in a Quantum Design PPMS-9 device. All 
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specific heat data measured in field were field cooled. 



III. RESULTS 

Fig.l shows the temperature dependent specific heat of 
annealed and unannealed Ba(Fei_ 2 -Coa;)2As2. As may be 
seen from the figure the annealing process has a dramatic 
impact on the specific heat of the samples measured. 
It increases T c between unannealed and annealed sam- 
ples quite substantially from 5.6 to 8 K (undcrdoped), 
20 to 25 K (optimaldoped), II to 17.2 K (overdoped, 
x = 0.105) and non-superconducting to 7.45 K (over- 
doped x = 0.14), where T c was determined by equal area 
construction of the specific heat jump. Moreover, the 
annealing process strongly suppresses the residual spe- 
cific heat 70 from 10.5, 3.6 and 14.6 mJ/mol K 2 to 1.3, 
1.3 and 3.8 mJ/mol K 2 respectively for under-, optimal- 
and over-doped samples. In the inset of Fig. lb we have 
also included data for an annealed optimally doped sam- 
ple from a different batch (sample s#2). This sample 
has T c = 26 K and the residual specific heat as low as 
0.25 mJ/mol K 2 which is the lowest reported value within 
the Co-doped Bal22 system. Furthermore, no Schottky- 
like anomaly is observed for the annealed materials. In 
addition to the increase of T c , the transition width broad- 
ens significantly for the overdoped sample, while it re- 
mains the same for the optimally-doped sample. For the 
underdoped annealed sample a reasonably sharp jump 
at 8 K exists in addition to a broad tail which starts at 
11 K, the temperature at which the resistive transition 
goes to zero (see Fig. 6). A similar broadening is observed 
in susceptibility measurements of Co-doped SrFe2As2^ 
and may indicate that the annealing process is not uni- 
form throughout the sample in some instances. 

As can be seen from Fig.l, above T c the C(T) curves 
overlap indicating very similar temperature dependence 
of the phonon contributions in all these materials. To 
obtain the electronic part of the specific heat, we have 
used the same approach used previously for the opti- 
mal and doped Ba(Fei_ x Co :r )2As2 samples (see RefH, 
HoUM). We have assumed that the phonon contribu- 
tion of the specific heat is independent of doping and 
we use the phonon specific heat obtained from the par- 
ent compound. Then we separate the lattice contri- 
bution to the specific heat of the pure BaFc2As2 com- 
pound as C ph (T) = C(T) BaFe2As2 
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7 BaFe 2 As 2T where 



^BaFe 2 As 2 j g the T _^ Q intercept of q j T of B a Fe 2 As 2 . 
So, the electronic specific heat is simply determined 

by C , e; (7 1 ) So (- Fe i-xCo JC ) 2 As2 _ Q^ji'\Ba{Fe- L - x Co a .)2As2 _ 

C{T) ph . 

Fig. 2 displays the electronic part of the specific heat 
of unannealed and annealed Ba(Fei_ a; Co a ;)2As2. In both 
cases it is obtained by subtracting the phonon contri- 
bution, together with a small Schottky-term below 1 K 
(for the unannealed sample), and normalizing by j n . A 
pronounced jump associated with the superconducting 
transition can be observed. As has been shown previ- 
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FIG. 2: (Color online) The temperature dependence of the 
electronic fraction of Ba(Fei_a; Co^ As2 normalized by nor- 
mal state specific heat (circles: x = 0.045, triangles: x = 
0.105, squares: x = 0.08). Empty and full symbols represent 
unannealed and annealed data respectively. 
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FIG. 3: (Color online) The low temperature dependence of 
the specific heat of the unannealed and annealed samples 
measured at different magnetic fields for H||c. The data are 
presented in the form of C/T vs. T 2 (see text) 
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FIG. 4: (Color online) (a) The magnetic field dependence 
of the low temperature specific heat of Ba(Fei_ I Co a ;)2As2 
(circles: x = 0.045, triangles: x = 0.105, squares: x = 0.08). 
Empty and full symbols represent unannealed and annealed 
data respectively, (b) the same data as above presented as 
A / y(H)/('j n — 70) vs. H/H C 2- The doted, solid and dashed 
lines described the field dependencies of the low temperature 
specific heat according to clean s— , anisotropic s— and clean 
d-wave^ models respectively. 

ously, none of the specific heat data can be described by 
a single s or d-wave mode ffi 10 i 38 ' 41 . This is also true for 
other Fe-based superconductors (see for instance Reflr3- 
H3). Clearly there are additional low energy excitations 
not captured by a single gap and, at least, a two gap sce- 
nario is necessary to describe the experimental results in 
this systern£ i 10 ' 38 i 41 . Furthermore, several curves have a 
hump at intermediate temperatures (for instance the an- 
nealed overdoped sample at 9 K) which can be taken as 
evidence for multiband SC as in MgB^. In addition to 



the suppression of the residual linear and Schottky terms, 
the finite temperature dependence of C e i/T for under- 
doped and overdoped samples changes with annealing, 
while there is no significant change at optimal doping. 
This can be seen clearly in the insets of figure 1. In 
particular, the low temperature dependence of the over- 
doped sample changed from C e ;/T = jq+BT 1 - 7 for the 
unannealed crystal to C e ;/T = -fo+BT 1 for the annealed 
crystal. For the undcrdoped annealed sample the tem- 
perature range over which the behavior C/T = jq+BT 1 
is visible is limited due to the lower T c , but the change 
in the temperature dependence at finite temperature is 
readily apparent from figure 2. In none of the samples 
do we find evidence for the C ~ J 13 / 2 behavior expected 
for samples with accidental points nodes in 2D 4 £ perhaps 
due to the finite number of impurities present in even the 
cleanest sample measured. 

The temperature dependence of the specific heat of an- 
nealed and unannealed samples, taken at different mag- 
netic fields applied along the c-axis is shown in Fig. 3. The 
field dependence of the residual linear term plotted in 
fig. 4a is extracted from this data by extrapolating the low 
temperature specific heat to K (from the temperature 
above the Schottky anomaly). A-y{H) = 70 (H ) - 70 (0T) 
reflects the rate of quasiparticle production with mag- 
netic field. For the unannealed optimally doped and 
overdoped samples (Figs 3c and 3e), aside from a broad- 
ening of the low temperature upturn with increasing 
field, C/T shifts rigidly upward as has been reported 
previously 8 -!*^. The effect of annealing has remarkably 
different effects on the field dependence of the two dop- 
ings however. For optimal doping the effect of a mag- 
netic field on the annealed crystal continues to be a rigid 
shift of C/T with increasing magnetic field, although the 
magnitude of A'j(H) is reduced by a factor of two rela- 
tive to the unannealed crystal. For the overdoped sam- 
ple however, the change in the residual linear term as a 
function of magnetic field A n f(H) grows by 40% in the 
annealed crystal compared with the unannealed sample. 
Furthermore, as we noted above the overdoped annealed 
sample has a temperature dependence of C/T= r jo+BT 1 
in the low temperature limit. The application of a 
magnetic field restores the temperature dependence to 
C/T=j (H) + BT 2 . The behavior on the underdoped 
crystals is qualitatively similar to that for the overdoped 
crystals. Again annealing enhances Aj(H) from 3.7 to 
8.8 mJ/mol K 2 in the undcrdoped sample. 

In an effort to determine whether the opposite evo- 
lution in the field dependence as a function of anneal- 
ing was due to a simple change in H C 2 or [7^-70 (0T)] 
we plot in figure 4b the field dependence of the resid- 
ual linear term of the specific heat scaled by [7,1-70 (0T)] 
versus H/H C 2- H C 2 was obtained from the specific heat 
determined T C (H) and applying the WHH formula 4 ^ to 
obtain H C 2(0), and is in reasonable agreement with the 
measured values of H C 2^ where available. For under 
and over-doping the field dependence of 7 scales with the 
unannealed crystals. However, at optimal doping there is 
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a dramatic departure from this scaling. The held depen- 
dence resembles the expectation for an isotropic s-wave 
superconductor, where the number of Caroli-deGcnncs 
bound states increases linearly with field due to the linear 
increase in the number of vortices entering the sample^. 
Alternatively, it would require an upper critical field of 
125 T to scale the data from the optimally doped an- 
nealed crystal with the remaining curves. This behavior 
is even more prominent in the annealed optimally doped 
sample with 7 =0.25 mJ/mol K 2 (sample s#2) where 
the data perfectly follows the prediction for an isotropic 
s-wave superconductor (see the filled stars in Fig. 4b). 
Consequently there exists a qualitative difference as to 
the effects of annealing between optimally doped and 
non-optimally doped crystals. A table summarizing the 
features described above are shown in Table I. 



IV. DISCUSSION 

The improved crystal quality of the annealed crystals 
as evidenced by the increase in T c , reduction in the resid- 
ual linear term of the specific heat, and suppression of 
any low temperature Schottky upturn puts us in a better 
position to try to understand the order parameter sym- 
metry of Co-doped BaFe2As2- We first discuss several as- 
pects of the data on a qualitative level, and subsequently 
we compare self-consistent calculations using a general 
two-band model which capture many of the salient fea- 
tures of the temperature dependent specific heat data at 
optimal and ovcrdoping. 



A. Gap phenomenology 

For the annealed optimally doped sample, where the 
residual linear term is the smallest, the zero field spe- 
cific heat of the annealed sample has a long temperature 
limiting behavior of C e i/T=j +BT 2 (see inset of Fig. 1). 
The value of B is 25 % larger than the T 2 term com- 
ing from the acoustic phonon branch. It is very difficult 
to get this precise low temperature dependence from a 
simple gap structure. Clean point nodes are capable of 
producing a specific heat C ~ T 3 , however, any impuri- 
ties which would be needed to provide the finite 70 term 
should also lift the point node and consequently remove 
the expected T 3 dependence. The fit which spans nearly 
a decade in temperature argues against this simply be- 
ing a crossover regime, although we can not rule it out 
entirely. Finally, bosonic modes might easily produce a 
C ~ T 3 contribution to the specific heat. Perhaps due to 
the suppression of the magnetic and/or structural tran- 
sition one may find such bosonic excitations which were 
not present in or modified from the parent 3 ^ or heavily 
overdoped compounds^. Regardless of the unknown ori- 
gin of the T 3 term at optimal doping, it is clear that in 
approaching the clean limit there are very few low energy 
electronic excitations. This is supported not only by the 



smallest residual linear term and the small increase in 
C/T with increasing temperature, but also by the fact 
that A7(9T) is the smallest in the optimally doped an- 
nealed sample. In fact, the normalized field dependence 
shown in figure 4b would indicate a nearly fully gapped 
superconductor. 

At optimal doping, prior to annealing, T c was lower 
indicating the existence of more impurities. In the 
unannealed sample C/T has a larger residual 7, and 
there also exists a larger field dependence to 7 relative 
to the annealed sample. These statements are consistent 
with a decrease of low energy excitations with annealing. 
A pure phase separation scenario can not explain both 
of these observations. Rather a gap structure which 
creates impurity states at low energies may be able to 
reproduce such an effect. 



The situation changes dramatically when moving to 
either under or over doping. We discuss the situation 
in going from optimal doping to ovcrdoping, and note 
that everything is qualitatively similar in going to under- 
doping as well. The annealed overdoped sample relative 
to the annealed optimally doped sample has several dif- 
ferences. The temperature dependence of the electronic 
specific heat for the overdoped samples is much better ex- 
pressed as C/T=7o+ST 1 (as opposed to C/T= l0 +BT 2 ) 
in which the magnitude of both terms is larger than that 
found at optimal doping. In addition, Aj(H) is larger, 
and is clearly sublinear. Furthermore, the temperature 
dependence evolves from C/T = 70 + BT 1 at zero field 
to C/T = jo+BT 2 at high fields. All of these features 
of the annealed over-doped sample are consistent with a 
gap structure that either has nodes, or very deep min- 
ima. Similar conclusions were previously found form 
a high quality overdoped crystals grown by a different 
technique 3 -^. 

The similarities and differences in the evolution of spe- 
cific heat features with annealing for the overdoped (un- 
dcrdoped) and the optimally doped materials is also in- 
structive (see tabic I). With annealing T c increases, 70 
drops, and a low temperature Schottky term disappears 
in all materials. On the other hand, in both overdoped 
and underdoped cases, the BT a and Aj(H) terms are 
enhanced in cleaner (annealed) samples, indicating that 
there are clearly more electronic excitations associated 
with the gap structure in the annealed case. This is pre- 
cisely opposite to our conclusion at optimal doping. The 
situation for the underdoped crystal is qualitatively iden- 
tical to the overdoped crystal. We cannot rule out the 
possibility that further improvement in sample quality 
for the overdoped and underdoped samples may lead to 
temperature and field dependencies which resemble the 
optimally doped sample. However, the qualitatively dif- 
ferent trends with annealing suggests that the gap struc- 
ture must qualitatively change as a function of doping. 
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TABLE I: Parameters obtained by fitting the equation C e i/T = 70 + BT a to the experimental data (below 5 K) of under 
(x=0.045), optimal (x=0.08) and overdoped (x=0.105 and 0.14) Ba(Fei_a;Coa;)2As2. Arrows mark the change (increase or 
decrease) of the parameter for the annealed sample compared to the unanealed one. 



sample [x] 


Tc [K] 


70 [mJ/molK 2 ] 


Schottky 


B [mJ/molK a+2 ] 


Q 


A 7 (9T) 


f)4J~i nn atitippiIpH 

\J . U LI 1 1 CXI 1 1 1 cL 1 L \_L 


5.6 


10.2 


No 


0.97 


1.2 


3.7 


0.045 annealed 


8.0 t 


1.0 4 


No 


3.71 t 


0.98 4 


8.8 t 


0.08 unannealed 


20 


3.6 


Yes 


0.12 


2.1 


6.07 


0.08 annealed 


25 t 


1.3 I 


No | 


0.09 4. 


2.2 t 


3.7 4 


0.08 annealed s#2 


26 t 


0.25 I 


No | 


0.04 4 


2.4 t 


3.85 4 


0.105 unannealed 


11 


14.5 


Yes 


0.47 


1.7 


6.1 


0.105 annealed 


17.2 t 


3.8 4. 


No | 


1.37 t 


1.05 4 


9.3 t 


0.14 unannealed 






No 








0.14 annealed 


7.45 t 


0.4 


No 


6.5 


0.9 


14.2 



B. Comparison to theory 

Currently the most highly discussed model for the su- 
perconducting gap in the iron pnictides is the so called s ± 
state which receives theoretical support from spin fluctu- 
ation models of various tight binding parameterizations 
of the Fermi surfaces of these materials^—. 

We computed the heat capacity within a model of equal 
size electron and hole Fermi surfaces with an interband 
pairing that supports an isotropic gap on the hole FS 
sheet, Ah(<p) = Ah and an anisotropic gap of the gen- 
eral form A e (<fi) = A e (l — r + r cos 2(j>) on the electron 
sheet, with parameter r varying from r = (isotropic 
state) to r — 1 (equal-size positive and negative lobes). 
The ratio of A e to A^ is determined by our choice of 
equal density of states on the hole and electron Fermi sur- 
faces. We accounted for the impurity scattering in the 
self-consistent T-matrix approximation, and solved the 
resulting gap equation to obtain the density of states and 
the entropy as a function of temperature. The impurity 
scattering is qualified by 7 = where T = J£™£ f) 

and rii mp is the impurity concentration. 

It is clear from the data that the optimally doped sam- 
ple exhibits significantly smaller temperature variation 
than either over- or under-doped materials, suggesting a 
very small density of low-energy excitations. Consistent 
with that observation, we find the best fit for these data is 
using a fully isotropic gap on the electron Fermi surface, 
r = 0, with relatively few impurities of moderate strength 
given by the phase shift 5 = 60 — 80° but significant in- 
terband scattering 5V = v^/va = 0.8 — 0.9, Fig.[5jt. For 
r = 0.0, S = 60°, 7 = 0.05 and 5V = 0.8, the calculation 
gives A e = -Ah = 3.8 meV. On the overdoped side, we 
found that several possible parameter values give an ad- 
equate fit to the data at the intermediate temperatures. 
However, the low-temperature behavior is most consis- 
tent with a nodal order parameter on the electron Fermi 
surface sheet (r > 0.5). For fully gapped systems, even 
when the impurity concentration is sufficient to produce 
a non-vanishing density of states at zero energy, we find 
significantly more features in the calculated C/T at low 
and intermediate temperatures, due to the strong varia- 



tion of the DOS between the upper edge of the impurity 
band and the superconducting gap edge^s. In contrast, 
in nodal systems the evolution of the Sommerfeld coeffi- 
cient is smooth as shown in Fig. [SJd, and the choices of 
r = 0.55 or 0.6, with closely spaced nodes on the electron 
FS sheets, give the best agreement with experiment. For 
r = 0.6, S = 60°, 7 = 0.15 and 8V = 0.6, the calcula- 
tion gives A e = 4 meV and Ah = -2.4 meV. The linear 
behavior down to the lowest temperatures suggests, from 
our fits, scatterers of intermediate strength, phase shift 
of close to 60°, and a substantial ratio of inter- to intra- 
band scattering. This is plausible in Co-doped systems 
where the dopants are on the Fe sites and directly af- 
fect the (i-orbitals contributing to both FS sheets. Using 
the same phase shift, but increasing the impurity con- 
centration and slightly increasing the relative strength 
of the interband scattering, we are able to fit well the 
unannealed data as well, as shown in Fig. \5]p. We also 
roughly account for the 64 % suppression of T c (17.2 to 
11 K) by increasing the impurity concentration from 0.1 
to 0.25 using r = 0.6, 5 = 60° and SV = 0.6 gives a 
62 % suppression of T c (from 0.77 T c0 to 0.48 T c0 ). The 
underdoped sample is in the region of the coexistence of 
magnetism and superconductivity. Consequently, we did 
not attempt a fit since our theory does not include the 
effects of the magnetic ordering that is likely significant 
in this regime. 



C. Effects of annealing 

Annealing clearly improves sample quality, as T c in- 
creases, the residual linear term decreases, and the low 
temperature Schottky anomaly is suppressed; however, 
the microscopic origin for these effects is less clear. The 
change of T c 's could suggest a change of cobalt concen- 
tration in the samples after annealing resulting in a shift 
in doping towards higher T c on the phase diagram. How- 
ever, this could be ruled out by the fact that T c increases 
for both the under-doped and overdoped samples. More- 
over, in the underdoped material, the spin density wave 
and structural transitions also increased (see Fig. 6) which 
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FIG. 5: (Color online) (a) Optimal doping, x = 0.08. Theo- 
retical fits (lines) to experimental data for annealed sample, 
(b) Overdoped regime, x = 0.105. Theoretical fits (lines) 
to experimental data for unannealed (open symbols) and an- 
nealed (filled symbols) sample. 



would contradict an increase in doping to account for the 
increase in T c . The increase in the spin density wave 
transition with annealing is also observed in the parent 
compound^ and in Co-doped SrFe2As2^. 

The presence of a finite residual linear term could orig- 
inate from a sample which contains both superconduct- 
ing and non-superconducting portions of a sample. Con- 
sequently, one might imagine that annealing is simply 
trading non-superconducting portions for superconduct- 
ing portions. However, as shown in Fig. 7, the volume 
of the superconducting fraction, as measured by mag- 
netic susceptibility, does not change significatively for 
the x = 0.105 sample. The fact that the overall mag- 
netic field dependence at optimal doping is suppressed in 
the annealed sample further disproves such a suggestion^. 



FIG. 6: (Color online) The temperature dependencies of nor- 
malized electrical resistivity of unannealed (squares) and an- 
nealed (circles) Ba(Feo.95sCoo.o45)2As2. Arrows mark the an- 
tiferromagnetic phase transition. Inset: low temperature part 
of the resistivity. 
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FIG. 7: (Color online) The temperature dependencies of 
magnetic susceptibility of unannealed (empty symbols) and 
annealed (filled symbols) Ba(Fen.895Coo.io5)2As2. The ZFC 
measurements have been performed on the same sample with 
the filed parallel to the ab plane of the crystal. 



Some fraction of the residual gamma may still originate 
from non-superconducting portions of the sample, but 
we believe the change as a function of annealing is most 
likely accounted for by a difference in the number of pair 
breaking impurity scatterers. However, we should cau- 
tion that annealing clearly does more than just change 
the number of pair breaking impurity sites, as evidenced 



8 



by its effect on the structural and SDW transitions. Also, 
the low temperature Schottky term is likely a result of 
free paramagnetic Fe or Co moments within the sample. 
Its suppression as a function of annealing indicates that 
these sites have either been removed from the sample, or 
possibly incorporated into the lattice. It is worth noting 
that a similar Schottky-typc anomaly was also observed 
in the cuprates but, up to now, the orgin of this effect 
remains unknown. Furthermore, the normal state Som- 
merfeld coefficient also changes with annealing from 13.7, 
18, and 23.2 mJ/mol K 2 prior to annealing to 14, 22, and 
20 mJ/mol K 2 after annealing for the under-, optimally-, 
and over-doped samples, respectively. 

While the temperature and field dependencies for the 
annealed samples now allow us to identify some differ- 
ence in low energy gap structure as a function of doping, 
it is interesting that the temperature and field depen- 
dence of all dopings with significant impurity concentra- 
tions could be scaled (after subtraction of the residual 
linear term) with one another—. This scaling persists in 
the temperature dependence of the annealed optimally 
doped sample^, but not its field dependence. However, 
the temperature dependence of the annealed under and 
overdoped samples no longer scale with the previous sets 
of data. 

All the results obtained clearly indicate that more work 
is needed to understand the microscopic effects of anneal- 
ing, and whether further improvement in sample quality 
can help elucidate additional details of the pairing mech- 
anism. 



V. SUMMARY AND CONCLUSIONS 

In summary, we report on the effect of annealing on 
the specific heat of under-, optimally- and over-doped 
Ba(Fci_ x Co x )2As2. We have shown that annealing at 
800°C for two weeks improves superconducting charac- 
teristics in the Co-doped BaFe2As2 system significantly. 
The heat treatment results in an increase of the critical 
temperatures and largely decreases the residual specific 
heat in all samples. Also, no Schottky-like term is ob- 
served at low temperatures. 

From the temperature and magnetic field dependence 
of the electronic specific heat of the annealed crystals 



we find a significant reduction of thermally excited low 
energy excitations for optimal doping relative to the un- 
der and over-doped regimes. This suggests a change in 
the gap structure as a function of doping which is con- 
sistent with the lifting of accidental nodes at optimal 
doping. Similar behavior has also been seen in Zn dop- 
ing studies^, thermal transport 2 ^, Raman scattering^, 
^tSRii and penetration depth measurements^. Whether 
the current theoretical models can self consistently ex- 
plain all these experimental observations remains to be 
seen. 

We note that the quality of samples has a strong in- 
fluence on the power law observed in the low tempera- 
ture specific heat data. Furthermore, the low temper- 
ature Schottky anomaly is also sample dependent and 
will likely manifest itself differently for different tech- 
niques. Consequently, the conclusions drawn by differ- 
ent techniques, and even identical techniques on differ- 
ent samples can dramatically vary. The power laws in 
transport and thermodynamic measurements are heav- 
ily relied upon to determine plausible gap structures and 
subsequently pairing mechanisms. It is also very difficult 
to have a single independent measure (such as RRR) to 
determine the quality of a sample. Consequently, this 
work emphasizes the need for multiple measurements to 
be carried out on the same high quality sample where 
possible. 
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